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Abstract 

This paper investigates the cuspidal configurations of 3-RPR parallel manipulators that may 
appear on their singular surfaces in the joint space. Cusp points play an important role in the 
kinematic behavior of parallel manipulators since they make possible a non-singular change 
of assembly mode. In previous works, the cusp points were calculated in sections of the joint 
space by solving a 24th-degree polynomial without any proof that this polynomial was the only 
one that gives all solutions. The purpose of this study is to propose a rigorous methodology 
to determine the cusp points of 3-RPR manipulators and to certify that all cusp points are 
found. This methodology uses the notion of discriminant varieties and resorts to Grobner 
bases for the solutions of systems of equations. 

KEY WORDS : Kincmatics, Singularities, Cusp, Parallel manipulator. Symbolic compu- 
tation 



1 Introduction 

Because at a singularity a parallel manipulator loses its stiffness, it is of primary importance 
to be able to characterize these special configurations. This is, however, a very challenging 
task for a general parallel manipulator. Planar parallel manipulators have received a lot of 
attention [1-5] because of their relative simplicity with respect to their spatial counterparts. 
Moreover, studying the former may help understand the latter. Planar manipulators with 
three extensible leg-rods, referred to as 3-RPR manipulators, have often been studied. Such 
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manipulators may have up to six assembly modes and their direct kinematics can be written 
in a polynomial of degree six [1, 2]. Moreover, they may have singularities (configurations 
where two direct kinematic solutions coincide). It was first pointed out that to move from 
one assembly mode to another, the manipulator should cross a singularity [3, 4]. Later, 
Innocenti and Parenti-Castelli [5] showed, using numerical experiments, that this statement 
is not true in general. In fact, this statement is only true under some special geometric 
conditions, such as similar base and mobile platforms [6, 7]. More recently. Macho et al. [8] 
proposed a method to plan non-singular assembly-mode changing trajectories. McAree [6] 
pointed out that for a 3-RPR parallel manipulator, if a point with triple direct kinematic 
solutions exists in the joint space, then the nonsingular change of assembly mode is possible. 
This result holds under some assumptions on the topology of the singularities [9]. For other 
mechanisms than 3-RPR manipulator, it is also interesting to note that encircling a cusp 
point is not the only way to execute a non-singular change of assembly mode [10]. A condition 
for three direct kinematic solutions to coincide was established in [6]. This condition was 
then exploited in [11] to derive a univariate polynomial of degree 96. A factored expression 
was obtained, one of the factors being a 24th-degree polynomial. The authors observed on 
many examples that the cusp points were each time defined by the 24th-degree polynomial, 
the remaining factors always defining spurious solutions only. However, they did not attempt 
to certify that the 24th-degree polynomial was really the only valid factor. Moreover, they 
never found more than 8 cusp points. 

The purpose of this paper is to propose a rigorous methodology to determine the cusp 
points of any 3-RPR manipulator and to certify that all cusp points are found. This method- 
ology uses the notion of discriminant varieties and resorts to Grobner bases to solve the 
systems of equations. Using symbolic computation, it is verified that the cusp points are 
really defined by a 24th-degree polynomial. For any given 3-RPR manipulator geometry, the 
maximum number of cusp points in sectional sections of the joint space is determined and 
the results are certified. In particular, a robot with 10 cusp points in a cross section of its 
joint space is found for the first time. 

The following section introduces the manipulators studied and recalls the main known 
results. Section 3 describes the algebraic tools used. Last section presents the methodology 
that is proposed to determine the cusp points using the algebraic tools. 

2 Modeling and state of the art 
2.1 Robot studied and Modeling 

A general 3-RPR planar parallel manipulator is shown in Figure 1. This manipulator has 
three extensible leg-rods actuated with prismatic joints. This manipulator can move its 
moving platform B1B2B3 in the plane. The vector of the three leg-rod lengths is L = 
{Pii P21 Pz) (Figure 1). The geometric parameters of the manipulators are the three sides of 
the moving platform di, d2, and the position of the base revolute joint centers Ai, A2 
and A^. The reference frame is centered at Ai and the x-axis passes through A2. Thus, 
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Figure 1: The 3-RPR parallel manipulator under study. 

Ai = (0,0), A2 = (v42^,0) and A3 = (^3^,^13^). 





Parameters 


Position variables 


Equation constraints 


MD 


P1,P2,P3 


01, 02, 03 


1 


\\B^B2\\ = d^ 
-82-83 = d2 

[ II 53^1 II =^3 


GRS 


Pl,P2,p3 


Bix, Biy, a 


1 


Pl^lll =Pl 
11^2^2 II = P2 

, Pa^all = P3 



Table 1: Two models for the 3-RPR manipulator 

The position of the moving platform can be expressed by two different sets of variables 
(see also Table 1): 

• MD Model (McAree-Daniel): the 3 angles 01,02,03- 

• GSR Model (Gosselin-Sefriou- Richard): the position coordinates of point Bi {Bi^, Biy) 
and the platform orientation a, = cos{a), ay = sin{a)). 

Each of these sets of variables leads to a different modeling system. 
2.1.1 MD Model 

The first model was used by MacAree and Daniel in [6] and has been extensively studied 
in [11]. Let 6 = (^1,6*2,^3) define the three angles between the leg-rods and the x-axis. 
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The six parameters {L, 6) define a configuration of the manipulator but only three of them 
are independent, so that the configuration space is a 3-dimensional manifold embedded in a 
6-dimensional space. The dependency between (i, 0) is obtained by writing the constraint 
equations of the manipulator, namely, the fixed distances between the three vertices of the 
mobile platform Bi, B2, B^: 

( T,{L, 0) = [b2(i, 0) - bi(L, 0)nh2{L, 0) - bi(X, 0)] - d\ 

r I r2(L, 0) = [b3(L, 0) - b2(L, 0)Y[h^{L, 0) - b2(i, 0)] - dl (1) 

[r3(L, 0) = [bi(L, 0) - b3(L, 0)Y\h^{L, 0) - b3(L, 0)] - dl 

where hi is the vector defining the coordinates of Bi in the reference frame as function of L 
and 0. 



2.1.2 GSR Model 

The GRS model was introduced in [12]. This model, combined with appropriate algebraic 
methods (see Section 3), will allow us to use a full symbolic computation of all cuspidal 
configurations. In this modeling, point Bi and angle a are used to specify the pose and the 
orientation of the platform. Let X = [Bi^, Biy,ax,ay), where (resp. ay) denotes cos{a) 
(resp. sin{a)). Let Px (resp. denote cos(/3) (resp. sin{P)). With these variables, we can 
parametrize the positions of the three points Bi, B2 and B^ as follows: 



bi 



Bly 



Bix + dittx 
Bly + diOy 



Bix + dsia^/Sx - ay/3y) 
Bly + d3{ax(3y + ay(3x) 



The lengths Pi, P2, Pa define the geometric constraints of the robot, and can be written 
as the euclidian norm of the vectors AiBi, A2B2, A^B^. Finally, the dependency between 
(-L, X) can be identified by writing the distances between the vertices of the moving platform 
and the vertices of the base platform. These distances pi, p2, Pa are the norms of the vectors 

A^i,A^2,A^3: 




= [bi(X, X) - ai]^[bi(Z, X) - ai] - p2 
= [b2(i, X) - a2]^[b2(i, X) - a2] - pI (2) 
= [h,iL, X) - a3]^[b3(i, X) - ag] - pi 

where a^ is the vector defining the constant coordinates of Ai in the reference frame. 
The formulas in Equations 2 can be expanded as follows: 

B,l + B,l-p\ = 

J. I (Bix + diax - A2xf + [Bly + dxoiyf - p\ =0 

^ {Bix + dsaxPx - d^ay/Sy - As^f + {Biy + d^a^Py + dsay/S^ - AsyY - pi = 

al + Oy-l =0 
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The last polynomial in E comes from the variable a replaced by the 2 variables a^,, ay 
and the equation constraint + — 1 = 0. This change of variables allows us to avoid the 
sine and cosine functions and keep the system algebraic. Moreover, this change of variables 
does not introduce any spurious solutions since the function 

[-7r,7r[ 

9 ^ (cos(^),sin(^)) 

is bijective on the circle C. 

The following subsections recall the main properties of the 3- RPR manipulator observed 
using the MD Model. 

2.2 Properties observed 

In a special configuration {L, 6), the manipulator meets a singularity, where two assembly 
modes coalesce. This happens when the constraint Jacobian drops rank. A 3-RPR parallel 
manipulator is in a singularity whenever the three lines {AiBi),i = 1..3, are concurrent or 
parallel [13]. 

When three assembly modes coalesce, the manipulator is said to be in a cuspidal confi- 
guration. Using series expansion, McAree and Daniel show that in such configurations, the 
manipulator loses first and second order constraints. This gives a necessary condition for a 
manipulator to be in a cuspidal configuration (see Table 4). 

Based on series expansion, [11] provided an algorithm to plot the singular curves in 
sectional slices of the joint space and to calculate the cusp points that appear on these 
curves. They calculated the cusp points in various sectional slices of the joint space for the 
manipulator defined by the following geometric parameters: 

Ai = (0, 0) A2 = (15.91, 0) A3 = (0, 10) di = 17.04 ^2 = 16.54 ds = 20.84 (4) 

The following facts were observed in [11]: 

• no more than 8 cusp points were found in a sectional slice defined for a given value of 
Pi; this result was also observed in many other manipulator examples; 

• the number of cusp points per section stabilizes to four when pi exceeds a certain value; 

• the set of all cusps points defined a set of curves in the 3-dimensional joint space 
(Fig. 2). 

2.3 Limitations of the previous study 

In [11], the cusp points were calculated in any sectional slice of the joint space by resorting 
to complex algebraic calculus based on expansion and elimination procedures. A univariate 
polynomial of degree 96 was obtained, which was shown to factor into several polynomials. 
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(a) A2x = 11, A3x = 7, A3y = 10 
di = 5, = 5, f3 = 37 
9^ = 53, 02 = 127, 03 = -90 
Pi = 5, p2 = 5, p3 = 3 
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Figure 3: A platform triangle 



(b) A2x = 11, A3x = 7, A3y = 10 
di = 5, (is = 5, /3 = —37 

= 53, ^2 = 127, 03 = -90 
Pi = 5, p2 = 5, p3 = 9 
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(left) and its mirror image (right) 



one of which is of degree 24 and was conjectured to be the one that gives all cusp points in 
the sectional slice at hand. The other factors were found to give always spurious solutions. 
The following remarks can be drawn: 

• there is no proof that the above mentioned 24-degree polynomial is indeed the only 
polynomial that gives all solutions; 

• the sectional slices were analyzed by a discrete method, and, thus, the analysis was 
not exhaustive; 

• a maximum number of 8 cusps points were found in many manipulator examples, but 
there is no proof that no more than 8 solutions can be obtained. 

• previous equation modeling does not distinguish between a platform triangle and its 
mirror image (see Figure 3), i.e. a triangle with the same di but with an angle /3 of 
opposite sign. 

One limitation of the previous studies came from the algebraic tools used to describe 
the robot, which introduced spurious components that had to be removed with empirical 
observations. In the literature, another method exists based on the elimination of variables, 
reducing the problem to the computation of the triple roots of an univariate polynomial 
[14, 15, 2]. Unfortunately, due to the size of the univariate polynomial, this method ran out 
of memory when applied to 3- RPR manipulators [11]. 
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The following section presents the minimal Discriminant Variety and the Jacobian crite- 
rion. These algebraic tools allowed us to compute and describe exactly the cuspidal config- 
urations of the 3-RPR manipulators, without introducing any spurious components. 

3 Algebraic tools 

This section recalls some mathematic definitions and their implementation in Maple soft- 
ware. The parametric systems we consider are supposed to have a finite number of complex 
solutions for almost all parameter values. This property is checked by computing a Grobner 
basis of the system at hand (see [16, page 274] and [17, Theorem 2] for more details). 

3.1 Discriminant Variety 

A Discriminant Variety V [17] is associated to a parametric polynomial system S. One of its 
main properties is that the real roots of S can be parametrized continuously on each open 
connected set in the complement of V. 



8 



Notations 1 

In this section, we denote by S{ti, ...,ts) a polynomial system of the form: 

Pi = 0, ...,Pm = and gi > 0, > 

depending on the s parameters ti, ...ts and the n variables Xi, ...,Xn- 

For more simplicity, (ti, ...,ts) will be omitted in the subsequent sections. 

For example, if we want to describe the number of solutions to the direct kinematics 
problem according to the articular parameters, then: 

• the polynomial equations pi, ...,Pm. are equations (3) with the specifications (4) 

• the polynomial inequations gi > 0, > are pi > 0, p2 > 0, > 

• the parameters ti are the articular variables pi,p2,P3 and the variables Xj are the pose 
variables Bi^, Biy,ax,(yy. 

The following definition is adapted from [17] for our problem. 
Definition 1 (Discriminant Variety) 

Let S be any parametric system. With the notation defined above, we call Discriminant 
Variety of S a set V of parameter values such that on each connected part ofW \ V: 

i. the roots of S can be parametrized as continuous functions of the parameters; 

a. the roots of S do not cross. 

Moreover, the intersection of all the discriminant varieties is called the minimal discrim- 
inant variety and is denoted by DV . 

Remark 1 The conditions i. and ii. imply that for every open set in the complement of 
DV , the number of solutions is constant. 

One of the main property of the minimal discriminant variety DV is that if S has finitely 
solutions for almost all parameter values, then for each connected cell U of W \ DV, the 
number of solutions of 5* restricted to U is constant. 

In the problem at hand, the minimal discriminant variety may be split into two compo- 
nents: 

- the inequation bounds in the parameters space. As far as a 3-RPR manipulator is 
considered, these bounds are trivial since the inequations already depend solely on the 
parameters; 

- the projection of the parallel singularities (defined in [18], see also section 3.3). These 
values are computed as follows: 
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1. Compute Ml, . . . , Mj-fc-j as the n x n minors of the Jacobian matrix: 

dxi dxn 



dpk dpk 
\ dxi dxn J 

2. Ehminate the variables Xi, . . . ,Xn in the system 

Pi = 0,. . = 0,det(M)i = 0, . . .,det(M)(^fc^) = 

The zero iso-surface of the returned polynomials is the projection of the desired 
singularities. 

3.2 Degree of a 0-dimensional system and multiplicity of a root 

The interested reader may find more detailed definitions on the degree and multiplicity in 
[19] or in [20, chapter 4]. 

Definition 2 (degree) Given a system S of polynomial equations with finitely many complex 
solutions, the degree of S is exactly the number of roots counted with multiplicities. 

The following unusual definition is exactly the same as the classical one [20] (dimension 
of a local ring), the reader will find the correspondence in [21]. 

Definition 3 (multiplicity) Let S be a system of polynomial equations in Xi,...,Xn with 
finitely many complex solutions. Let / : C" — )■ C 6e a linear form and P{Z) a univariate 
polynomial returned by the computation of a rational univariate representation [21]. 
Then P can be written as: 

p{z)= n iz-i{c^)r- 

a root of S 

We call multiplicity of a the exponent fia • 
Remark 2 When S is reduced to a univariate polynomial equation in X with k roots 

k 

1=1 

then the multiplicity of each root Ci is exactly /x^. 
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By extension, we can also define the multiplicity of the roots of a parametric system. 

Definition 4 (multiplicity in a parametric system) Let S{ti, ...,ts) be a parametric system 
depending on the parameters ti, ...,ts and on the variables Xi...,Xn- Assume furthermore that 
for all parameter values {t^, G C^, the number of complex roots of S{ti, is finite. 

Then, the multiplicity of a root a = {t^, ...,t^,Xi, ...,xD is the multiplicity of the root 
{x^, in the system S{t1, 

Remark 3 This multiplicity definition can be naturally extended to define the multiplicity 
of a leaf above a connected component in the complement of the discriminant variety. 

3.3 Jacobian criterion 

We recall here the Jacobian criterion, used in this article to extract the points of multiplicity 
greater than or equal to 2. More details can be found in [22]. The Jacobian matrix to 
be considered here is exactly the same as the one used in [18] to define the singularities of 
type 2 of a manipulator. In this article, we will refer to these type-2 singularities as parallel 
singularities. 

Theorem 1 Let pi, . . . ,pm be parametric polynomials of C[ti, ...,ts][xi, . . . , Xn] admitting 
a finite set of common complex roots for each value of the parameters ti, ...,ts. Let J be the 
Jacobian matrix related to the variables Xi, 

dpi dpi 
dxi dxn 



dpm dpm 

dxi dxn 

// Ml, . . . , M^m^ denotes the n x n minors of J, then for all t° = (t^,...,t^) G C^, the 
following parametric system: 

Pl(t°,Xi, ...,Xn) = 0, . . .,Pm{t°,Xi, ...,Xn) = 

det(M)i(t°,xi,...,Xn) = 0,...,det(M)^^^)(t°,xi,...,x„) = 

has a solution in Xi, ...,Xn if and only if the system 

Pl(t°,Xi, ...,Xn) = 0, . . . ,Pm{t^,Xi, ...,Xn) = 

has at least one root of multiplicity greater than or equal to 2. 

This theorem can be used to define singular configurations as in [18], but also to define 
the cuspidal configurations, as shown in section 4.2. 
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3.4 Software implementation in Maple 



The algebraic tools used in this paper were adapted from the SALSA library developed at 
INRIA. The manipulations and computations of the algebraic objects are performed with 
the math software MAPLE. The functions used are described below. 

FGb Given a set of polynomials generating an ideal /, f gbrs : -f gb_gbasis allows the 
computation of a system of generators of the ideal. These generators provide a way to 
reduce any polynomial to a canonical unique form modulo J. With the option "elim"=true, 
it makes it possible to compute an algebraic system, the zero of which is the closure of the 
projection of the roots of /. 

This elimination of variables is based on Grobner basis computations, using the algo- 
rithms F4 ([23]). As explained in the following sections, this function is used to compute the 
singular and cuspidal positions, as well as the discriminant variety. 

RS A system S of polynomial equations pi = 0, . . . ,pk = is said 0-dimensional if its 
number of complex solutions is finite (this can be straightforwardly tested on a Grobner 
basis of the system). 

In this case, the RootFinding: -Isolate function computes directly the real solutions of 
S. By default, each solution is given by a box with rational bounds, ensuring that no pairs of 
boxes overlaps. The algorithms implemented in Maple are based on the Rational Univariate 
Representation ([21]). 

This function allows us to compute the number of cuspidal positions. 

Discriminant Variety Given a parametric system of equations pi = 0, . . . ,pk = and 
inequations gi ^ 0, . . . , Qr ^ 0, ti, . . . ,ts being the parameters and Xi, . . . ,Xn the unknowns, 
its minimal discriminant variety is computed with the algorithm presented in [17]. 

As shown further in section 4.3, the discriminant variety is a key tool for the description 
of the cuspidal positions. 

This function is part of the package RootFinding [Parametric] in MAPLE. 

Limitations These functions stand only for polynomials with rational coefficients. In 
particular, the trigonometric functions sin, cos, tan have to be replaced by algebraic variables. 

4 New modelling and method 

In this section and the subsequent ones, we use the equation system and the notations 
induced by the GSR Model presented in section 2.1.2. In this model, the position of the 
platform is given by 4 variables: Bi,^,Biy the coordinates of the point Bi and ax,(yy, the 
cosine and sine of the angle a respectively. 
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The configurations of a 3-RPR manipulator may be classified into three categories: the 
regular configurations, the singular configurations and the cuspidal configurations. The de- 
scription of the special configurations (singular and cuspidal) is important for path planning. 
First, the regular configurations are the most common ones: these are the configurations 
where the position of the platform is locally uniquely determined by the values of the control 
parameters. Then the singular (resp. cuspidal) configurations are two (resp. three) coa- 
lesced configurations. These configurations were studied through local analysis in the past, 
but they can be also described with algebraic methods, as shown further in this article. 

4.1 Singular configurations 

From a geometrical point of view, we saw in Section 2.2 that a singular configuration repre- 
sents coalesced assembly modes. 

From an algebraic point of view, this condition can be retrieved in term of multiplicities. 
When the variables L is specialized with real numbers, the set of equations (3) has finitely 
many solutions. Moreover, a multiplicity is associated with each solution (see Section 3.2). 
Using this notion, it can be verified that the singular assembly modes are exactly configu- 
rations that are solutions of multiplicity greater than or equal to 2. Indeed, the Jacobian 
criterion [22] states directly that these configurations are parallel singularities. 



4.1.1 Modelling 

To describe the singular points, we use the notion of multiplicity in a 0-dimensional system. 
The system of equations (3) is 0-dimensional when the lengths pi,p2,p3 are given: it has a 
finite number of complex solutions. In this case, each solution can be associated with an 
integer greater than or equal to 1: its multiplicity, as defined in section 3.2. 



Definition 5 A configuration Pq = {p'l, p2, p'^, Bi^, Biy, a'^, ay) of the 3-RPR manipulator 
is said singuL 
equal to 2 of: 



is said singular if and only if (i?i|^., a^, a^J) is a solution of multiplicity greater than or 



E((p;,p°,p°),X) = [ ] 



T 



To find the singular configurations, we use the Jacobian criterion, reminded in section 3. 
Applying the Jacobian criterion [22] on the polynomial system Equations (3), we deduce 
that the singular configurations are exactly the solutions of the system: 



EfL,X) 







T 



(5) 



Remark 4 By adding the determinant of the Jacobian matrix to Equations (3), we get a 
new system whose solutions are included in the original system. More precisely, the roots of 
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Equations (5) are exactly the roots of multiplicity greater or equal to 2 of Equations (3). The 
solutions to Equations (3) of multiplicity 2 are solutions of multiplicity 1 to Equations (5), 
and the solutions of multiplicity greater or equal to 3 in Equations (3) have a multiplicity 
greater or equal to 2 in Equations (5). 

Moreover, using Grobner basis computations to eliminate variables with Maple , a poly- 
nomial in the parametric variables L, denoted Psingi can be computed. Its solutions are 
exactly the sets of length values for which the manipulator admits singular configurations. 

4.1.2 Example 

In this section, the computation results are given for the manipulator example presented in 
[5, 6, 11, 24]. 

Once substituting the numerical values of Equations (4) into Equations (5), a system is 
obtained, the solutions of which are exactly the singular configurations of the studied ma- 
nipulator. After eliminating the variables X with Grobner bases computation, a polynomial 
Psing in L is obtained of total degree 24. The zeros of this polynomial define a surface in 
the parameter space. By specifying pi = 14.98 as in [11, 6], the curve as in [11, 6] can be 
observed in Figure 4. 

4.2 Cuspidal configurations 

Cuspidal configurations are associated with second-order degeneracies that appear for triply 
coalesced configurations. As shown in [6, 11], these configurations play an important role in 
path planning. 

To find cuspidal configurations, the idea of [6, 11] was to analyze the kernels of the 
matrices in the first and second order terms of the series expansion of Equations (1). It 
allowed the authors to find the cuspidal configurations automatically when pi is specified 
as a numerical value. However the description they get is a list of slices of the cuspidal 
configurations, and they miss what could happen between 2 successive slices. Using the 
notion of discriminant variety and a generalization of the Jacobian criterion, we introduce 
a complete certified description of the cuspidal configurations. In particular, our approach 
allows us to find cuspidal configurations that were missed in the previous papers, and to 
certify that all cuspidal configuration are determined. 

4.2.1 Modelling 

Like for the singular points, we define the cuspidal configurations algebraically. 

Definition 6 The configuration Pq = {pi, p2, P3, -Bi°, of the 3-RPR manipulator 

is said cuspidal if and only if a^., dy) is a solution of multiplicity 3 to 

E((p;,P^,P^),X)=[ ]^ 
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In Section 4.1, the Jacobian criterion allowed us to select the configurations of multi- 
plicity higher than or equal to 2. However, we now want to compute the configurations of 
multiplicity 3. Fortunately, we saw in Remark 4 that these configurations were the roots of 
multiplicity 2 of Equations (5). We can thus use the Jacobian criterion on Equations (5) to 

get a system whose roots of multiplicity 1 describe exactly the cuspidal configurations. 

dB 

Let det(J) = det(— — ) and let Mi, M2, M3, M4, M5 be the 4x4 minors of the Jacobian 
O'-X. 

matrix of the 5 polynomials Ei, E2, E3, E4, det(J) with respect to the variables X. Then 
the following system of 9 equations defines exactly the cuspidal configurations of the 3- RPR 
manipulator: 



(6) 



E(L,X) = [00 
det(J) = det(M)i 

det(M)2 = det 



det(M)3 = det 



det(M)4 = det 



det(M)5 = det 



Remark 5 These equations are algebraically dependent. In particular, even if the number 
of equations is 9 and the number of variables is 7, the set of solutions is a curve of dimension 
1 in 






V 






= 








dEi 


dE2 


dE, 


9 det (J) 


9X 


dX 


dX 


dX 


dEi 


dE2 


9E4 


(9 det (J) 


dX 


dX 


dX 


dX 


dEi 


dEs 


9E4 


(9 det (J) 


dX 


dX 


dX 


dX 


dE2 


dE3 


9E4 


9 det (J) 


dX 


dX 


dX 


dX 



System (6) contains equations of degree smaller than 5. This will allow us to compute a 
certified description of the cuspidal curves in Section 4.3. 



4.2.2 Example 

For the robot example defined by Equations (4), triple points are exactly the solutions of 
Equations (6). Moreover, when pi is specified. Equations (6) has finitely many solutions. 
Using the methods of real solving for 0-dimensional systems of the function Isolate of 
Maple, one can find easily the roots of Equations (6) for any given pi. 
The system defined by Equations (6) has 9 equations: 

• the first four equations are the geometric model equations, each being of degree 2; 

• the 5th equation is the Jacobian of these equations and has degree 3; 
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P2 


P3 




Bly 




ay 


1 


0.845 


3.777 


5.336 


-13.997 


0.633 


0.773 


2 


13.851 


6.260 


-14.963 


0.698 


0.998 


-0.045 


3 


31.276 


16.178 


-6.104 


13.679 


-0.543 


-0.839 


4 


17.988 


26.446 


14.721 


-2.769 


-0.985 


0.167 


5 


30.449 


26.619 


-10.363 


10.816 


.537 


0.843 


6 


16.027 


29.566 


14.437 


3.995 


0.999 


-0.010 



Table 2: The 6 cuspidal configurations (roots of Equations 6) for pi = 14.98 



• the last four equations are obtained through the iteration of the Jacobian computation, 
each one is of degree 5.. 

The full equations can be found at the address http://www.irccyn.ec-nantes.fr/ 
~chablat/3RPR_cuspidal .html. 

In particular, for pi = 14.98, we get 6 triple points (see Table 2), which confirms the 
results of [11]. 

Moreover, in [11], the authors have shown that for different numerical values of pi, among 
different spurious factors, it is possible to extract a univariate polynomial in ti := tan{6i) of 
degree24. 

Using our cusp modeling, if we add the new variable ti = tan(6'i) = Biy/Bi^ and elimi- 
nate the variables p2, ps, Bi^, BiyOx, ay from Equations 6 with Grobner basis computation, 
we can compute (in 9 seconds) a bivariate polynomial Q in pi and ti. The real roots defined 
by Q is the closure of the projection of the cuspidal configurations. We can see that Q is 
irreducible and has a degree 24 in ti. In particular, by substituting a numerical value of pi 
in Q, we get directly a univariate polynomial of degree 24 in ti. This generalizes and certifies 
the results in [11]. 

4.3 Cusps analysis 

Equations (6) define implicitely the cuspidal configurations of the robot. Without further 
computations, these equations are not sufficient to describe the geometry of the triple roots 
of the 3- RPR manipulator. In [6, 11], the authors observed that the set of cuspidal con- 
figurations is finite in slices of the robot joint space for given values of pi. This leads to 
the conjecture that the set of cuspidal configurations forms a curve in the full joint space 
(Pi! P2i Ps)- Furthermore, the authors of [11] described the number of cuspidal configurations 
for discrete values of pi in M. This allowed them to find configurations that were omitted in 
previous works. 

In this section, we give a certified and exhaustive description of the number of cuspidal 
configurations as function of pi. In particular, this study allows us to find values of pi, for 
which the 3-RPR manipulator has 10 cuspidal configurations, while previous works never 
observed more than 8 cuspidal configurations. This gives new possibilities to change assembly 
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mode without crossing singularities. More generally, our method compute the number of 
cuspidal configurations for all the values of pi in open intervals, and not only slices. 

The dimension of the set of solutions of Equations 6 in the complex field is 1. This can 
be verified by computing and analyzing the Grobner bases of Equations 6 (see [25, chapter 
9] for more details). 

To describe geometrically the solutions of Equations 6 as function of pi, we consider it 
parametric system where: 

• the single parameter is pi 

• the unknowns are p2, ps, -Bi^^, a^, oty 

First, using [16, page[341], we can check that for almost all values of pi, the roots of 
Equations 6 have multiplicity 1 and are thus corresponding exacty to cuspidal configurations. 

Then we follow the work of [17, 15, 26] to describe the roots of a parametric system. Our 
process has 2 steps: 

• we first compute its minimal discriminant variety. 

• we then compute the number of solutions for sample parameter's values chosen outside 
the discriminant variety 

4.3.1 Discriminant variety of the cuspidal configurations 

We consider Equations 6 as a system parametrized with pi. In this case, its minimal dis- 
criminant variety is a finite set of values of pi denoted by: 

DV = {ai,...,afc} C M,A; > 

The main property of the discriminant variety is that for each value of pi in an interval 
]aj,aj+i[, the Equations 6 has the same number of distinct real solutions. 

The points of the discriminant varieties are the real roots of univariate polynomials. 
The lines pi of table 3 show numerical approximations of these roots. The full univariate 
polynomials are not given here for lack of space. They can be found at the address http : 
//www. irccyn. ec-nantes . f r/~chablat/3RPR_cuspidal . html. They were computed in 11 
seconds with a 2.9GHz Intel cpu. 

The lines ^Cusp gives the number of cuspidal configurations for pi in an interval ]aj, aj+i[, 
where the are the real values of pi defining the discriminant variety DV. 

4.3.2 Number of cuspidal configurations 

Using the main property of the discriminant variety, we know that to count the number of 
cuspidal configurations outside the discriminant variety, it is sufficient to compute a finite set 
of sampling points in the complement of DV and to solve the corresponding zero-dimensional 
systems. More precisely, we choose a value Vi in each interval: 

joj, aj+i[, such that Oj, Oj+i are consecutive reals of DV 
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#Cusp 
Pi 


2 4 2 
0.000] [ 0.148 ] [1.655] [ 1.660 ] [ 


#Cusp 
Pi 


4 6 8 6 
2.261] [ 2.975 ] [9.186] [9.186*] [ 


7^ Cusp 
Pi 


8 6 8 6 
9.257 ] [ 9.257* ] [10.905] [10.905*] [ 


#Cusp 
Pi 


8 6 8 6 
14.579] [14.579*] [20.555] [ 20.562 ] [ 


#Cusp 
Pi 


8 10 8 6 
26.786] [ 28.094 ] [28.107] [ 28.257 ] [ 


#Cusp 
Pi 


8 6 4 


30.740] [ 30.779 ] [30.946] 



* these values are close but distinct, the closest values differ by at least 10 

Table 3: Discriminant variety of the cuspidal configurations (Equations 6) w.r.t. pi; the 
numerical values were truncated 



19 



and compute the number of solutions to Equations 6 when pi = Vi. Moreover, we choose a 
value Vk+i in ja^, +oo[, and count the number of solutions of Equations 6 when pi = Vk+i- 
The results of these computations are summarized in the lines T^Cusp of Table 3. To count 
the number of real solutions of Equations 6 when pi = Vi, we use the real solver given by 
RootFinding: -Isolate in Maple. It took around 5 minutes to solve the systems induced 
by each of the 23 sample points with a 2.9GHz Intel cpu. 

From Table 3 we conclude that, asymptotically, the manipulator has 4 cuspidal configu- 
rations, and that this number is stable as soon as pi is greater than 31. Moreover, we can 
observe that the robot may have up to 10 cuspidal configurations when pi is in the interval 
]28.095, 28.107[. Figure 5 shows the singular curve of the parallel robot for pi = 28.10 and 
the corresponding 10 cuspidal configurations. 

5 Conclusion 

This paper shows that efficient algebraic tools can be applied to analyse in a certified way 
important kinematic features of parallel manipulators such as the determination of cusp 
points. These points are known to play an important role in planning non-singular assembly 
mode changing motions. 

A new method was introduced, which is able to characterize all the cusp points for the 
3-RPR manipulators. It allowed us to determine the number of cusp points for all the slices 
of the joint space. For the first time, a section for a given pi was found with 10 cusp points. 

Table 4 summarizes the method used in this article and relates it to the 2 previous main 
methods to compute cuspidal configurations in the literature. 

This results can be computed with a toolbox available now in Maple 12 and later versions 
[27]. 
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Figure 5: Singular curve for pi = 28.10. The circles show the 10 cuspidal configurations in 
this configuration. 
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